function plotPositiveDomain(model,elements)
ne = model.countElements();
mesh = model.getMesh();
for k=1:ne
    if elements(k).isCrack
        PosDomain = elements(k).PosDomain;
        np = size(PosDomain(:),1);
        f = mesh.getFace(k-1);
        nodes = f.getNodes();
        for i=1:np
            X = nodes(PosDomain(i)).getPosition();
            text(X(1),X(2),'+','color','red');
        end
    end
end

end
